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Abstract 

A new approach based on antisymmetrized molecular dynamics is proposed to correctly take 
account of the many-body correlation. We applied it to the spectra of low-lying, light and strange 
baryons. The inclusion of the quark-quark correlation is vital to predict the precise spectra, and 
the semi-relativistic kinematics is also important to correct the level ordering. The baryon spectra 
calculated by the present method is as precise as the Faddeev calculation. 
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Although recent lattice calculations based on Quantum Chromodynamics (QCD) start 
providing reliable hadron spectra, it is still difficult to obtain precise, detailed predictions 
of the physical states [ll . Thus, QCD- inspired, effective models are still useful to get some 
insight into many phenomena of hadronic systems. The central issue to be addressed is then 
the quantitative description of low-energy phenomena, e.g., hadron spectra, baryon-baryon 
interactions, in- medium change of hadron properties [2|, etc. 

Among QCD- inspired models, the simplest approach is the nonrelativistic (NR), con- 
stituent quark model (CQM), and it is doubtlessly successful in describing the ground-state 
properties of hadrons and baryon-baryon interactions at low energies . There are, how- 
ever, some obscure problems such as the level-ordering problem in excited hadron spectra, 
etc. In addition to ordinary mesons and baryons, the (NR) CQM is also applied to exotic 
states, e.g., multiquark states including pentaquarks, non gg-mesons (tetraquarks 3]), etc. 

The QCD lagrangian for massless up and down quarks is chirally symmetric, and the axial 
symmetry is spontaneously broken as evidenced in the absence of parity doublets in the low- 
mass hadron spectra. This implies the existence of the massless Nambu-Goldstone (NG) 
bosons, e.g., the pions. The non-zero pion mass is then a consequence of the fact that the 
light quark has a small mass, which gives the explicit symmetry breaking. Thus, one arrives 
at a low-energy scenario that consists of the NG bosons and constituent quarks (with a 
rather massive, constituent quark mass in the NG phase) interacting via a force governed by 
spontaneously broken, approximate chiral symmetry. Thus, it is of very importance to take 
into account the effect of the NG-boson exchange as well as the usual one-gluon-exchange 
(OGE) force in the CQM. Note that multigluon degrees of freedom can be elminated by 
introducing a confining potential. Furthermore, the relativistic effect may be vital to produce 
the hadron spectra, because the quarks move in the region of the hadron size. 

There are a lot of model calculations to study hadron spectra, baryon-baryon interactions, 
etc jsl, y, y, 0, [s]. In the calculation of baryon spectra, the most precise method may be 
based on the Faddeev equations for the three constituent quarks j3] . However, for multiquark 
states like the exotic hadrons, it is quite difficult to perform such precise calculations. Thus, 
it is very important and useful to construct a powerful method which allows us to precisely 
calculate the wavefunction even for multiquark states. 

Molecular dynamics (MD) is very successful in treating many-body systems. In particular, 
fermionic molecular dynamics (FMD) was ffist developed by Feldmeier jol to describe the 
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ground states of atomic nuclei and heavy ion reactions in the energy regime below particle 



production. Antisymmetrized molecular dynamics (AMD) 10|, is very similar to FMD 



with respect to the choice of the trial state. In the past decade, AMD has been applied 



to various studies of light or medium nuclei 



12| . In AMD, it is not necessary to take any 



model assumption (like shell, cluster, etc.), and the AMD wavefunction can simultaneously 
describe a variety of nuclear structure by following the variational principle. 

However, AMD has not yet been successfully applied to a few-body system, because the 
trial wavefunction in AMD is given by the Slater determinant of single-particle, gaussian 
wave packets and hence the correlation among nucleons is missing. To improve this weak 
point, some of the present authors have recently proposed a modified version of AMD, in 
which the Jacobi coordinates and the generator-coordinate method (GCM) are introduced 



to describe the correlation 



13j.^ We here refer to this method as the Jacob i-coordinate-basis 



AMD (JAMD). In JAMD, one can easily treat the many-body correlations. Furthermore, 
it is possible to extract the center-of-mass (cm.) wavefunction and remove the zero-point 
energy. 

In this paper, we apply the present method to the SU(3) octet- and decuplet-baryon 
spectra, and compare the JAMD result with that of the simple AMD or the Faddeev result 
to illustrate how this method is useful. 

In AMD, the wavefunction of A-quark system with definite (±) parity, |\&^), is given 



by Qm 



(1) 



i=l 



where l^^p^), \^Pi) and \C) are, respectively, the spatial, spin, flavor and color wave- 

functions. The color wavefunction is then antisymmetrized as 



\c) = ^sgn(p,)|m;;^(i)m;;^ 



(2) 



■ m. 



i=l 



where m'^ specifies the color and sgn(pj) is the sign of the permutation pi {i = 1, 2, 



Pi 



1 2 3 ... A 



Jl J2 33 ■■■ 3 A 



1 2 ^ ... A 

p,(l) p,(2) p,(3) ... p,{A) 



(2) 



A!): 



(3) 



^ The nuclcon-nuclcon correlations are also considered in the stochastic variational method (SVM) 
or the coupled-rearrangement-channel Gaussian-basis variational method (CRCGV) 15 1. 
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Because the color wavefunction is already antisymmetrized, the other wavefunctions must 
be symmetrized in Eq.([T]). The spatial wavefunction is then given by the product of a 
single-particle wavefunction 

= ii±p)\MMm2ipm---MMm. (4) 

where 1 ± P is the parity projection operator and the i-th single-particle wavefunction (at 
Tj) is given by a gaussian function 

= exp[-iy{fi - Zj)\ (5) 

with Z the center of the wave packet and v a variational parameter for its width. The spin 
wavefunction is expressed as 



E 



%^m■■■p^{A)\^p,(l) 



■ m: 



(6) 



where the coefficient, c^, is a variational parameter. However, since, for the low-lying baryon 
states, the spin-flavor structure (|i5p.)|jFpJ in Eq.([T])) can be given by SU(6) symmetry, we 
here use it to reduce the computation time. 

In JAMD, the spatial wavefunction is expressed in terms of gaussian functions with the 

n 

Jacobi coordinates i = 1,2, ■ ■ ■ , A) pj[. Then, the usual single-particle coordinates, rj, 
can be related to the Jacobi coordinates through pi = X]fc=i -^jfc^fc; where the matrix K is 
given by [l^l 

/ 1 -1 ■■■ \ 



K 



nil 
mi2 



m2 
mi2 



m2 



-1 







V 



•"112- -A-i rni2...A-i 
n m2 



-1 



(7) 



— / 

■A ' 



where vnyi-.-k = f^ii + "^2 + ■ ■ ■ + fnk and is the mass of the i-th particle. The center of 
the gaussian wave packet, Xi, in the JAMD wavefunction is also related to Zj using the K 



matrix 



13|. Thus, the spacial part of the JAMD wavefunction, (^■^ , with the width, v^,, can 



be expressed in terms of the usual coordinates as 



exp 



exp 



A A 

[ - K^jiPi - ^i) ■ {Pj - ^j) = exp [ - ^ Uijifi - Zi) ■ {fj - Zj) 



(8) 



where each element of the width matrix, Uij, is treated as a variational parameter, and ufj, 



Uij and Uij = 1,2, ■ ■ ■ , A) satisfy the following relation: 



k,l=l 



k=l 



(9) 



Now let us apply the JAMD method to the SU(3) octet and decuplet baryons. The 
hamiltonian is given by 



H = T + V^{fij) + VoGEirij) + VcoN (fij). 



(10) 



where T is the kinetic-energy term. In this paper, as well as the NR form, we consider the 
semi-relativistic (SR) form to take into account the relativistic effect: 



1 sr^A rrii+mj 



{Pj - Pif + Y^t=i n^i for NR, 



for SR, 



(11) 



where rrii is the i-th quark mass. 

The potential generated by the exchanges of the NG bosons (vr, K, rj) and the a meson is 
given by V^x(^»i) ~ ^7r(^«j) + Va{rij) + Vxirij) + , where 

A2 |- A3 -, 3 

^ y (m^r,,) - -^YiA^r,,)\ (a, ■ a,) ^(A," ■ AJ), (12) 

a=l 



9ch 1^1 



4:71 12m,mi — ml 

J TV TV 



Vain,) 



4tt Ki — ml 



A. 



Y{m^rij) -Y{A„r 



(13) 



9ch 



A2 



4:71 Vlmimj K\ — m\- 



Y{mKn,) - ^Y{AKnj) (a.-a,)^(A,"-A»), (14) 



a=4 



9!h 



A? 



Att 12mimj A^ — 



mi 



(15) 



with Tij = Vi — fj, A'* the SU(3) generator, gch the chiral coupling constant, A^^o-./r.r; the cutoff 
parameter, Y{mr) the Yukawa function and Op(= cos^^p(Af ■ A^) — sin^p) the mixing term 
for cous,deriug the physical „. The masses of the NG bosons are denoted by ,n, arrd 
m^. The a mass is chosen to be ~ + 4m^^ [T^. 



It is possible to treat the LS or tensor force in the JAMD method [3| • 
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The potential due to the OGE is [g, [3] 



VoGE{rij) = - asXi ■ Xj h - agXi ■ Xj 




1 1 r r 1 e-'^'^/'^o 



h for NR, 

/ \ (16) 

l[l + la,-a,j forSR, 

with the quark-gluon couphng constant and tq a parameter to be fixed from the data. 
Note that in Ref. [7| ro is chosen to be a function of the quark mass. However, in the SR 
calculation, we can obtain a good result of the baryon spectra even if ro is assumed to be 
constant (see Table HTl) . 



The confining potential is given by 



m 



acrjj(A--Ap for SR, 

with Oc the effective confinement strength, /ic a parameter for the screening effect and Csg 
a parameter for taking account of the vacuum energy. 

To perform the numerical calculation, we use the following identities: 

1 1'°° 1 

Jf^m^ = VC-^ / -^(e~^P'+'^'^^' -e-^'"')dw, (18) 

e~4;;^~^ (iw, (19) 





where C is an arbitrary, real number. Using these formulas, one can expand the SR kinetic 
energy and the potentials in terms of gaussian functions. 

Varying the variational parameters included in the JAMD wavefunction, the total energy 
of the system 

is minimized. To perform such calculation, it is very convenient to use the frictional cooling 
method which provides the time-development equations for the center, Zi, and the 



6 



width parameter, Uij, as^ 

^ = /^ff^ (z = l,2,...,A; A = x,y,z), (21) 

^""IS- ('>:' = 1.2.-,^). (22) 

where /i and /i' are arbitrary, negative real numbers. After sufficient time steps for the 
coohng, we can obtain the optimized wavefunction as a function of Zi and We refer to 
this wavefunction as v|/J^^i^-i. 

In contrast, we consider another wavefunction which is given by superposing the J AMD 
wavefunction ^ (± is suppressed here) with Zi calculated in vI/JA'^d-i ^^^^ with Vij supplied 
by a geometrical progression or a random number generator [l^. In this case, Zi and 



are no longer the variational parameters. Then, the wavefunction (we call this ^J^^^^ iij 
is given by 

^JAMD-II ^ ^^^^^^ ^ ^^^2,^ ^ ^ . . . ^ (23) 

where, instead of Zj and z/j^, the coefficients, c, c", c'", ■ ■ ■ , are now variational parameters, 
and they are determined by the Hill- Wheeler equation 

5((^JAMD~II|^|^JAMD-II^ _ ^^^JAMD-II|^JAMD-II^^) _ q_ (24) 

Now we are in a position to show the numerical result for the baryon spectra. The 
parameters in the present calculation are listed in Table [H In case of the NR calculation, 
we take the parameters given in Ref. T]. We assume that the value of tq in Voge is flavor- 
independent. In the SR calculation, referring to Ref. 6j, we determine the parameters in 
the potentials. Finally, three parameters, Cs^, as and m^, remain. Then, the nucleon (iV) 
mass is reproduced by tuning Csg, while as is chosen so as to fit the N-/S. mass difference. 
The strange-quark mass, m^, is determined from ffi to the N-VL mass difference. 

In Table [Tll we present the results of AMD and JAMD, and compare them with the 



Faddeev calculation 



7|. Furthermore, in Fig. [H the deviation of the calculated, ground- 



state mass from the observed one is shown to illustrate the model dependence of the baryon 



^ Because the Jacobi coordinates, {pi,Xi), can be transformed into the old set, {ri,Zi), using the K 
matrix, it is sufficient to know the time development of the latter variables [l3| . The actual calculation is 
performed under the constraint of Xa = 0. 
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TABLE I: Values of the parameters. 









NR [T] 




SR [6] 


Fixed 












u, d quark mass 






313 




313 


NG bosons 


TTi-Tr ffm^"'^) 

' ' "TV V / 




0.70 




0.70 




m„ (fm^^) 

'/ V / 




2.77 




2.77 




niK (fm~^) 




2.51 




2.51 


cr meson mass 






3.42 




3.42 


cutoff 


Ajr = Arr (fm~^) 




4.20 




2.20 




A„ = (fm~^) 




5.20 




2.70 


coupling constant 


5,\/(4vr) (fm-i) 




0.54 




0.54 


mixing angle 






-15 




-15 


confinement 


ar (MeV) 




230 




110 








0.70 






OGE 


' u 




0.35 




0.74 








NR 




SR 






AMD 


JAMD-I 


JAMD-II 


JAMD-I JAMD-II 


Free 












s quark mass 


rus (MeV) 


598 


587 


554 


562 525 


vacuum 


Cs, (MeV) 


333 


346 


372 


128 177 


OGE 




0.858 


0.761 


0.540 


0.775 0.500 



masses. The present result clearly show that the JAMD approach is much better than 
AMD in describing the baryon spectra. This fact implies the importance of the quark-quark 
correlation in the baryon structure. In particular, within the NR calculation, the JAMD-II 
result is very close to the spectra given by the Faddeev approach. 

In the SR calculation, the present approach again reproduces the spectra of the low-lying 
light and strange baryons very well. In particular, the JAMD-II can provides the precise 
result as in the Faddeev calculation. It is also remarkable that the level ordering of the lowest 
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TABLE II: Calculated baryon masses (in MeV). The excited baryons are calculated in the JAMD- 
II. 



state 






NR 




SR 


Experiment 




AMD 


JAMD-I 


JAMD-II 


Faddeev|Y] 


JAMD-I 


JAMD-II 




A r / -1 /n-i- \ 

N{l/2^) 


939 


939 


939 


939 


939 


939 


939 


N{l/2-) 


1481 


1460 


1409 


1411 


1522 


1480 


1553 


N {1/2^) 






1423 


1435 




1409 


1440 


A(3/2+) 


1230 


1232 


1234 


1232 


1235 


1236 


1232 


A*(3/2+) 






1589 






1601 


1600 


S(l/2+) 


1315 


1264 


1246 


1213 


1244 


1214 


1193 


S*(l/2+) 




1400 


1678 


1644 




1553 


1660 


S(3/2+) 


1435 


1404 


1396 


1398 


1384 


1382 


1385 


S(l/2-) 


1722 


1639 


1631 


1598 


1726 


1678 


1620 


A(l/2+) 


1166 


1160 


1135 


1122 


1139 


1120 


1116 


H(l/2+) 


1418 


1401 


1375 


1351 


1374 


1359 


1318 


1^(3/2+) 


1673 


1673 


1673 


1650 


1671 


1670 


1672 
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> 

J 80 
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2 60 
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20 
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FIG. 1: Difference between the calculated mass and the experimental value. 
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FIG. 2: Root-mean-square radius. 



positive- and negative-parity states in the nucleon spectra can be correctly reproduced in the 
SR calculation ja, ^. This fact certainly results from the relativistic kinematics, because, in 
the NR calculation, we cannot produce the correct level ordering. It is, however, necessary 
to study further to obtain a quantitative result of N{l/2~) and N*{l/2^). 

In Fig [21 we show the root-mean-square (rms) radius of the baryon calculated using the 
quark wavefunction. It is noticeable that the rms radius is very small in the SR calculation, 
whereas, in the NR calculation, it has the reasonable size. It should, however, be noticed that 
the meson cloud surrounding the baryon core contributes to the present value additionally. 
The tendency of the rms radii of the low-lying baryons does not depend much on the choice 
of the model. 

In summary, we have proposed a new approach of molecular dynamics based on AMD, in 
which the many-body correlation can be considered correctly. In this paper, we applied it to 
the spectra of the low-lying, light and strange baryons. It is shown that the inclusion of the 
quark-quark correlation is very vital to predict the precise spectra, and that the relativistic 
effect is also important to correct the level ordering. The JAMD approach can reproduce the 
baryon spectra very well. In particular, the spectra given by the JAMD-II calculation are 
very similar to the Faddeev result. Although we study the ordinary baryons in this paper, 
the present approach is very promising even in the calculation of a system containing four. 
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five or more quarks. Thus, it is very intriguing to calculate the spectra of exotic hadrons 
like pentaquark, tetraquark, etc |l7l |. 
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